† Corresponding author. E-mail:
To improve the modeling accuracy of radiative transfer, the scattering properties of aerosol particles with irregular shapes and inhomogeneous compositions should be simulated accurately. To this end, a light-scattering model for nonspherical particles is established based on the pseudo-spectral time domain (PSTD) technique. In this model, the perfectly matched layer with auxiliary differential equation (ADE-PML), an excellent absorption boundary condition (ABC) in the finite difference time domain generalized for the PSTD, and the weighted total field/scattered field (TF/SF) technique is employed to introduce the incident light into 3D computational domain. To improve computational efficiency, the model is further parallelized using the OpenMP technique. The modeling accuracy of the PSTD scheme is validated against Lorenz–Mie, Aden–Kerker, T-matrix theory and DDA for spheres, inhomogeneous particles and nonspherical particles, and the influence of the spatial resolution and thickness of ADE-PML on the modeling accuracy is discussed as well. Finally, the parallel computational efficiency of the model is also analyzed. The results show that an excellent agreement is achieved between the results of PSTD and well-tested scattering models, where the simulation errors of extinction efficiencies are generally smaller than 1%, indicating the high accuracy of our model. Despite its low spatial resolution, reliable modeling precision can still be achieved by using the PSTD technique, especially for large particles. To suppress the electromagnetic wave reflected by the absorption layers, a six-layer ADE-PML should be set in the computational domain at least.
To improve the precision of climate modeling and atmospheric remote sensing, radiative transfer models (RTMs) that can accurately calculate the radiation transferred through the atmosphere with aerosols and clouds are required.[1–3] As fundamental input parameters for the radiative transfer simulation, the light-scattering properties of aerosols (especially those in the solar spectrum) should be accurately modeled.[4–7] However, owing to the irregular shapes and inhomogeneous compositions of aerosol particles (such as mineral dust and soot), their light-scattering processes are not adequately understood, and substantial uncertainties still remain in their optical properties.[8,9] To simplify the scattering process of aerosol particles in the RTMs widely used at present, nonspherical aerosol particles are usually taken as spherical ones with equivalent volumes or surface areas, which will definitely decrease the computational accuracy of radiative transfer.[6,10–12] Many researchers have also found that nonspherical aerosol particles exert a significant influence on the polarized components of radiation.[6,9,13]
With the awareness of the importance of the scattering properties of nonspherical particles, many scattering calculation models have been established.[8,14,15] According to the size parameter ranges applicable to these models, they can be classified into three theoretical branches: the models in the Rayleigh regime, semi-classical scattering models, and models in the geometric-optics (or physical optics) regime. Models in the Rayleigh regime and geometric-optics regime include the Rayleigh approximation, the Rayleigh–Gans–Stevenson approximation (RGA), the anomalous diffraction theory, the bridging method,[16,17] and the geometric optics approximation (GOA).[8,18] Limited to their physical essence, they are mainly applicable to particles with sizes much smaller or larger than the light wavelength.[8] However, because the size of aerosol particles is comparable to the light wavelength in the solar spectrum, these models cannot be applied to calculate their light-scattering properties. For semi-classical scattering models, a particleʼs scattering properties are obtained by solving Maxwellʼs curl equations and Helmholtz equations analytically or numerically. Models of this type can be grouped into three categories, i.e., field expansion models, volume integral methods, and microelement models. The T-matrix method (TMM),[19,20] separation of variables method (SVM),[21] and point matching method (PMM)[22] belong to the first category, where TMM has become a widely tested model for the improvements made by Mishchenko, Bi, and other researchers.[8,19,23–25] By using these models, high calculation precision can be easily achieved with low memory consumption and high computational efficiency.[8] However, owing to the limitation of the boundary conditions for numerical calculation, most of these models can only be applied to particles with regular shapes. Moreover, their numerical simulation processes become nonconvergent for particles with large sizes, large refractive indices, or extreme shapes. The volume integral methods include the method of moments (MoMs)[26] and discrete dipole approximation (DDA).[27–29] The main advantage of DDA and the MoMs is that they can be applied to particles with any irregular shapes and inhomogeneous compositions. However, it should also be noted that although simulation accuracy will be improved with increasing dipoles, the matrix equations that need to be solved will become larger as well. Therefore, to ensure stability and efficiency of the computational process, DDA and the MoMs are mainly applicable to small particles.[30] With the development of computational electromagnetics, such micro-element models as the finite difference time domain (FDTD) and finite element method (FEM) have gradually been introduced into the light-scattering simulation of nonspherical aerosols.[31,32] The main advantage of these methods is that they can overcome the singularity of volume integral methods; thus, they are more suitable for particles with complex shapes and large sizes. However, owing to the constraints of numerical dispersion and the Courant stability condition, the spatial resolution of the computational domain should be smaller than λ/10 (where λ is the light wavelength); therefore, the time and memory consumed would be huge for the light-scattering simulation of large particles.[4]
To overcome this problem, the pseudo-spectral time domain (PSTD) was introduced into the numerical simulation of the light-scattering process by Liu and Yang.[33] In this model, the Fourier-transform-based pseudo-spectral method and the finite difference technique are used to calculate the spatial and temporal derivatives of Maxwellʼs equations, respectively.[34] Because of the excellent performance of the Fourier transform, this method can achieve very high simulation accuracy with a relatively coarse spatial resolution, and therefore it can decrease the computational time and memory dramatically.[33] Despite the excellent performance of PSTD, there are two aspects that need to be improved: first, owing to the specificity of its discretization form, the incident field is mainly introduced by a pure scattering field technique in the PSTD scheme now, which is difficult and time-consuming for particles with large size and complex shapes; second, the absorption boundary conditions (ABCs) employed in the PSTD scheme are mainly limited to Berengerʼs perfectly matched layer (BPML) and uniaxial perfectly matched layer (UPML), whose performance is relatively low compared to convolution perfectly matched layer (CPML) and ADE-PML for the FDTD scheme.[35] To improve modelʼs performance, Sun et al.[4] generalized the CPML for the PSTD model. Here, we generalize the weighted total field/scattering field (TF/SF) technique[36] and ADE-PML for the 3D-PSTD scheme and an improved PSTD scattering model of our own is established. In our scattering model, not only can the simulation errors caused by the wave reflection of ABC be reduced, but the incident wave can also be introduced without considering the actual shapes and sizes of the scatters.
The structure of our paper is as follows: in Section
The basic structure of the PSTD scattering model is presented in Fig.
![]() | Fig. 1. (color online) Basic structure of the parallelized PSTD light-scattering model for nonspherical aerosol particles. |
The preprocessing module is employed to prepare the input parameters for the light-scattering computation. These input parameters include the propagation direction and polarization states of incident light, the refractive index of the aerosol particle, and the particleʼs shapes discretized by a grid mesh. The electromagnetic field computational module is designed to calculate the electromagnetic field in the near- and far-field regions. There are three key techniques in this module, i.e., the ADE-PML technique for the truncation of the computational space, the weighted TF/SF technique for the introduction of the incident wave and the near-to-far transformation technique, where the near-to-far transformation is performed in the frequency domain, while the others are implemented in the time domain. The scattering parameter calculation module is applied to calculate a particleʼs scattering parameters based on the near and far electromagnetic field components; these scattering parameters include the absorption and extinction cross-sections, the scattering phase matrix and single scattering albedo.
In the PSTD model, the scattering process is simulated by numerically solving Maxwellʼs equations, where Maxwellʼs equations in an isotropic medium have the following form:
![]() |
![]() |
![](cpb_27_5_054215/cpb_27_5_054215_ieqn1.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn2.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn3.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn4.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn5.gif)
In contrast to the FDTD scheme, in the PSTD model, the pseudo-spectral method is applied to calculate the spatial derivatives of Maxwellʼs equations, whereas the time derivative is discretized by a central differential scheme. The basic principle of the pseudo-spectral method is to use Fourier integrals to obtain the derivatives of the field components, where the Fourier integrals are carried out by using the fast-Fourier-transform (FFT) algorithm. Because of the high-precision reconstruction ability of the FFT/IFFT technique, the spatial derivatives can be calculated accurately with a relatively coarse grid mesh. According to the idea of the PSTD scheme, the derivatives of the electromagnetic components along the x-, y-, and z-axes can be presented as
![]() |
![]() |
![]() |
![](cpb_27_5_054215/cpb_27_5_054215_ieqn6.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn7.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn8.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn9.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn10.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn11.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn12.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn13.gif)
Because of the abrupt changes in the optical properties at the surface of aerosol particles, the electric and magnetic field will be discontinuous at the interface between the scatterer and the ambient medium, which will result in unnecessary high-frequency components in the FFT process (i.e., the “Gibbs phenomenon”). In order to suppress the generation and accumulation of the unnecessary signals, the high-frequency components in the frequency spectrum are directly truncated before the IFFT process, which is similar to the method adopted by Liu et al.[33] Taking the derivatives along the x-axis for example, after the improvement, the calculation equation for
![]() |
![](cpb_27_5_054215/cpb_27_5_054215_ieqn15.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn16.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn17.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn18.gif)
Based on the theoretical analysis above, Maxwellʼs equations can be easily discretized by replacing the spatial derivatives with
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![](cpb_27_5_054215/cpb_27_5_054215_ieqn20.gif)
![]() |
![]() |
Although the electromagnetic scattering process is an open-space problem, the calculation of the PSTD should be always implemented in a finite computational domain owing to the limitations of the computer memory and CPU. Therefore, approximate ABCs have to be imposed to truncate the computational domain in order to suppress wave reflections back into the domain and permit all outward-propagating wave analogs to exit the domain, almost as if the domain were infinite.
In our model, the ADE-PML, an excellent ABC used in the FDTD scheme is generalized and applied to the 3D-PSTD scheme. The main advantage of the ADE-PML is that its implementation is completely independent of the ambient medium and high absorption performance can be achieved without increasing the memory consumption. Similar to CPML, the ADE-PML can be derived from modified Maxwellʼs equations in the frequency domain[37]
![]() |
![]() |
![](cpb_27_5_054215/cpb_27_5_054215_ieqn21.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn22.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn23.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn24.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn25.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn26.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn27.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn28.gif)
Substituting
![]() |
![]() |
Furthermore, transforming Eqs. (
![]() |
![]() |
Similar to the method adopted in Section
![]() |
![]() |
![]() |
![](cpb_27_5_054215/cpb_27_5_054215_ieqn35.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn36.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn37.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn38.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn39.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn40.gif)
![]() |
The implementations of the ADE-PML for the other electromagnetic field components can be similarly obtained, and thus their derivations are not presented here.
The introduction of an incident wave is another key technique in the PSTD scattering model. Because the spatial distribution of the electromagnetic field components is not defined by Yee cells, the traditional TF/SF technique for the FDTD scheme is no longer applicable to PSTD. Therefore, in our model, the weighted TF/SF technique proposed by Gao[36] is generalized and applied to the 3D-PSTD scheme.
Similar to the traditional TF/SF technique, the weighted TF/SF technique introduces an incident wave by adding incident terms into the connection region between the total field and the scattered field. As shown in Fig.
![]() |
In Eq. (
The wave propagation equations in the connection region can be derived by substituting
![]() |
The incident wave can be properly introduced by solving the equations above. However, from these equations, two problems remain to be dealt with in order to solve the equations above numerically, i.e., the selection of the window function ζ and the discretization of the modified Maxwellʼs equations.
In the weighted TF/SF technique, the selection of the window function is very important, because weighing the incident field with window functions in time domain means the convolution operation between the incident wave and window function in frequency domain. In our PSTD model, the IBH function (the integral form of the Blackman–Harris function) is set as the window function (the IBH function is widely applied to the digital signal processing field and shows the best performance[36]).
To meet the requirements of sampling accuracy, the number of sampling points in the connection layer should equal or larger than 2(M+1). For example, if M is set as 3, then an eight-layer connection region should be set at least. Because the computational space of the PSTD scattering model is three-dimensional, then the window functions should be expanded to 3D space correspondingly. In our model, the 3D window function is constructed by the product of the window functions in different directions, expressed as
![]() |
After the window function is determined, another problem is the discretization of the modified Maxwellʼs equations. Similar to the method introduced in Section
![]() |
![]() |
![](cpb_27_5_054215/cpb_27_5_054215_ieqn47.gif)
Similar to the FDTD scheme, the calculation of the scattering amplitude matrix and phase matrix should be based on the far electric fields obtained for two incident waves with orthogonal polarization states. Therefore, the surface integral method based on the Huygens principle is employed to transform the near-field components to their far-field in our model.
In the surface integral method, the far-field scattered by the aerosol particles is regarded as the field radiated by the equivalent electric and magnetic currents at a closed surface in the scattering field region with the scatterer inside. According to the electromagnetic equivalence theorem, the equivalent electric and magnetic currents can be calculated from the electric and magnetic components at the extrapolation surface
![]() |
![]() |
![](cpb_27_5_054215/cpb_27_5_054215_ieqn48.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn49.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn50.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn51.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn52.gif)
With the equivalent electric and magnetic currents obtained, the far electric field can be calculated by the following equations:
![]() |
![](cpb_27_5_054215/cpb_27_5_054215_ieqn53.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn54.gif)
![]() |
![]() |
![](cpb_27_5_054215/cpb_27_5_054215_ieqn55.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn56.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn57.gif)
![]() |
![](cpb_27_5_054215/cpb_27_5_054215_ieqn58.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn59.gif)
![]() |
![](cpb_27_5_054215/cpb_27_5_054215_ieqn60.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn61.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn62.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn63.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn64.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn65.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn66.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn67.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn68.gif)
The integral scattering properties of aerosol particles include absorption cross-sections, extinction cross-sections, scattering cross-sections and single scattering albedo. These parameters reflect the total ability of the aerosol to absorb and scatter electromagnetic waves. The calculation models for these parameters are similar to those designed for the FDTD and the MRTD,[14] where the absorption cross-section Cab can be computed by normalizing the energy depleted by aerosols with the total energy of the incident wave[14]
![]() |
![](cpb_27_5_054215/cpb_27_5_054215_ieqn69.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn70.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn71.gif)
The extinction cross-section Cext can be calculated by the volume integral of the extinction Poynting vector
![]() |
![](cpb_27_5_054215/cpb_27_5_054215_ieqn73.gif)
![](cpb_27_5_054215/cpb_27_5_054215_ieqn74.gif)
The PSTD scattering model for nonspherical aerosol particles was coded by Fortran 90, and the OpenMP technique was employed for the parallelization of the model. To quantitatively evaluate the performance of the PSTD scattering model, the results of this model are compared with those obtained from the well-tested models such as the Lorenz–Mie theory, Aden–Kerker theory and TMM in this section for spherical particles, inhomogeneous particles and nonspherical particles.
The scattering phase matrices simulated by the PSTD scattering model is first validated against Lorenz–Mie theory for two spherical particles, a small one and a large one.
Figure
As can be seen from the figures, the nonzero phase matrix elements obtained by PSTD closely agree with the exact solutions given by the Lorenz–Mie theory well, where in forward directions (
The precision of the PSTD model is further validated for a large spherical particle with a size parameter of nearly 100. In this test, the light wavelength and the refractive index are the same as those for the small one; the radius of the particle is set as
From the figure, it can be found that the results of the PSTD scheme show excellent consistency with those of the Lorenz–Mie theory in terms of their overall variation patterns, illustrating that our PSTD model can simulate the scattering properties of large particles with reasonable precision. For the scattering phase function, its simulation errors fluctuate with scattering angle notably with their values ranging from near-zero in forward directions to 38.5% at the scattering angle near 180°. For
To validate the PSTD modelʼs ability to calculate the integral scattering properties of aerosol particles, the extinction efficiency Qext, absorption efficiency Qab and single scattering albedo ω are calculated for particles with different sizes, and the results are listed in Table
![]() | Table 1.
Comparison of the integral scattering parameters simulated by Lorenz–Mie theory and the parallelized PSTD scattering model. . |
For the case of inhomogeneous particles, the PSTD model is validated against the Aden–Kerker theory for a sphere with a core–mantle structure. In this test, the wavelength of incident light is set as λ =
The curves of the phase matrix elements obtained by the PSTD model approximately coincide with those of the Aden–Kerker theory, where the relative simulation errors of F11 are less than 10% for most scattering angles, and the absolute errors of
The precision of the PSTD model is also investigated for nonspherical particles by comparing the scattering phase matrix and integral scattering parameters calculated by the PSTD model with the exact solutions of the T-matrix method.
To illustrate the performance of the PSTD in phase matrix calculation, the nonzero phase matrix elements of a cylindrical particle are given in Fig.
From the figure, it can be easily found that excellent agreement is achieved between the results of the PSTD model and the TMM, validating the precision and reliability of our model. For the scattering phase function, its simulation errors are generally smaller than 10% at scattering angles ranging from 0° to 90°; in the backward scattering directions, the simulation accuracy is reduced to some extent, but its simulation errors are smaller than 15% for most scattering angles. The angular distribution of the simulation errors of
To further test the modeling accuracy of the PSTD model for integral scattering properties, the extinction efficiencies (Qext), absorption efficiencies (Qab) and scattering efficiencies (Qsc) of the cylinders are also calculated by the PSTD model and the TMM, respectively, and the results are listed in Table
![]() | Table 2.
Comparison of the integral scattering parameters calculated by the PSTD model and the T-matrix method. . |
It can be found that the relative errors of the extinction and scattering are all smaller than 1% for particles with different shapes and refractive indices, and the errors of the calculated absorption efficiencies are generally smaller than 2%, indicating the high calculation accuracy of the PSTD model for nonspherical particles. On the whole, the extinction efficiencies obtained by the PSTD model are larger than those of the TMM, while the absorption efficiencies are generally underestimated for most cases. From the variation of the simulation errors, it can be found that the simulation accuracy of the integral scattering properties is gradually reduced with the increase of refractive index.
The model is further validated with the DDASCAT developed by Draine and Flatau.[28] In this process, the particle is set as a spheroid with a spherical core, where the horizontal axis and the rotational axis of the spheroid are set as
As can be seen from these figures, the phase matrices obtained by DDASCAT and the PSTD model show good consistency, where the scattering phase functions obtained by the two models almost coincide with each other and their relative discrepancies are less than 10% for all directions; for
The integral scattering properties obtained by the PSTD model are also compared with those obtained by DDASCAT for three inhomogeneous nonspherical particles, the results are presented in Table
![]() | Table 3.
Comparison of the integral scattering parameters calculated by the PSTD model and DDASCAT. . |
Spatial resolution is not only an important factor determining the construction accuracy of a particleʼs shape, but it is also a parameter influencing the discretization accuracy of Maxwellʼs equations. Owing to this reason, the influence of the grid size on calculation accuracy is systematically investigated in this section. To quantitatively evaluate the modeling precision, we use the following parameters to evaluate the performance of the PSTD model, i.e., the relative simulation errors of integral scattering parameters, the square root error of the scattering phase function and linear
![]() |
![]() |
Table
![]() | Table 4.
Influence of grid size on the simulation accuracy of aerosol scattering properties. . |
As shown in the table, with the increasing spatial resolution, the modeling accuracy of the phase function and integral scattering parameters are improved as well on the whole, where for particles with a size parameter of 10, in case that the grid size increases from λ/8 to λ/20, RSEF11 decreases from 11.6672% to 6.7534%, and RSEF12 is reduced by 0.1025. It can also be found that even if the spatial resolution is as low as λ/8, very high calculation accuracy can still be achieved by the PSTD model, especially for large particles. This phenomenon can be explained by the high-precision reconstruction ability of the pseudo-spectral method, because the FFT/IFFT technique is used in the PSTD model, very high precision can be achieved even with a spatial resolution near the Nyquist sampling limit, thus the spatial derivatives can be calculated accurately with a relatively coarse grid mesh. When the spatial resolution is given, the modeling accuracy is reduced as the increase of particle size. According to simulation errors for particles of different size, if we take the required precision criteria as follows: (i) the relative error of Qext is less than 1%, (ii) RSEF11 is less than 25%; then we can find that, for particles with
The performance of ABC is also an important factor influencing modeling accuracy, because the electromagnetic field reflected by ABC layers can alter the exact distribution of the field components. Generally, the thicker the ADE-PML is, the better the absorption ability of ADE-PML can be achieved, but the computational burden will be increased as well. Therefore, it is necessary to analyze the ADE-PML performance and determine the optimal thickness of the ADE-PML for the PSTD scattering calculations.
For this purpose, the simulation accuracy of a scattering phase matrix for ADE-PML layers with different thickness is first investigated and their performance is compared with that of a six-layer BPML (BPML is widely applied in the traditional PSTD model). In this simulation, the light wavelength is set as
As can be seen, with the increasing thickness of the ADE-PML, the calculation accuracy is improved correspondingly, where when the thickness of the ADE-PML is larger than six layers, the error curves are almost coincident with each other, indicating that the modeling accuracy is not affected by the electromagnetic reflection of the ADE-PML. Simulation errors were mainly caused by the starting approximation of particle shape and the discretization of Maxwellʼs equations. Therefore, we can conclude that an ADE-PML with a thickness larger than six layers is enough for a scattering simulation. Comparing the results of the six-layer ADE-PML and BPML, it can be found that although their variation patterns are similar, the simulation errors of ADE-PML is slightly smaller than that of BPML, indicating that the performance of ADE-PML is better than BPML in PSTD calculation.
The performance of the ADE-PML was also investigated for integral scattering properties, and the results are demonstrated in Table
![]() | Table 5.
Influence of the thickness of ADE-PML layer on the modeling accuracy of the integral scattering properties (spheroid case). . |
The performance of the ADE-PML is further evaluated for spheroidal particles, and the results are shown in Fig.
The parallel computational efficiency of the PSTD model is evaluated in this section. In this test, the light wavelength is
![]() | Table 6.
Parallel computational efficiency of the PSTD scattering model. . |
From the table, we find that the computational time of the PSTD model is reduced dramatically after the parallelization design, where for a particle with a size parameter of 20, the parallel acceleration ratio Sp reaches 3.8148 when six processors are used. It also can be found that η decreases as the number of processors increases; this can be explained by the increase of the communication time between different processors. Compared with the traditional PSTD model, it is also found that the improved PSTD model can reduce the computational time to some extent, where for the particle with a size parameter of 40, the time needed can decrease by approximately 9%.
The uncertainties of the scattering properties of aerosol particles with irregular shapes and inhomogeneous compositions are an important factor limiting the modeling accuracy of radiative transfer. To this end, a light-scattering model based on the PSTD technique is presented. Using this model, scattering by particles with arbitrary shapes can be effectively simulated. In this model, the ADE-PML, an excellent absorption boundary condition for FDTD, is generalized and applied to the PSTD scheme; the weighted TF/SF technique is employed to introduce an incident wave into the 3D computational domain; and to obtain the electromagnetic field in the far region, the surface-interpolation method based on the Huygens principle is used for the near-to-far transformation. The precision of the PSTD model is validated by comparing its results with those of the Lorenz–Mie theory, TMM, DDASCAT and Aden–Kerker theory, and the influence of the spatial resolution and thickness of the ADE-PML on the modeling accuracy is analyzed as well. Several conclusions are drawn as follows.
(i) The calculation results of the PSTD model, irrespective of the scattering phase matrix or the integral scattering parameters, show great consistency with those of the well-tested scattering models, indicating the high simulation accuracy of our model.
(ii) With the decreasing of the spatial resolution, the calculation accuracy of the PSTD models is improved as well. For particles with size parameters smaller than 40, a grid size of λ/24 is fine enough for light-scattering simulation, whereas for particles with size parameters larger than 40, a spatial resolution of λ/16 is enough.
(iii) The performance of the ADE-PML is closely dependent on its thickness. To suppress the wave reflections back into the domain, the thickness of the ADE-PML should be larger than six layers.
[1] | |
[2] | |
[3] | |
[4] | |
[5] | |
[6] | |
[7] | |
[8] | |
[9] | |
[10] | |
[11] | |
[12] | |
[13] | |
[14] | |
[15] | |
[16] | |
[17] | |
[18] | |
[19] | |
[20] | |
[21] | |
[22] | |
[23] | |
[24] | |
[25] | |
[26] | |
[27] | |
[28] | |
[29] | |
[30] | |
[31] | |
[32] | |
[33] | |
[34] | |
[35] | |
[36] | |
[37] | |
[38] |